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ABSTRACT 

This study on the Coma cluster suggests that there are deviations from a completely random galaxy 
orientation on small scales. Since we found a significant coincidence of hot-gas features identified 
in the latest X-ray observations of Coma with these local anisotropies, they may indicate regions of 
recent mutual interaction of member galaxies within subclusters which are currently falling in on the 
main cluster. 

Subject headings: galaxies: clusters: individual (A1656) — galaxies: evolution — galaxies: interactions 
— galaxies: statistics — X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy orientation studies in the past (see 
lHawlev fc Peebles! 119751 and references therein) were 
chiefly carried out in order to verify miscellaneous 
galaxy cluster formation scenarios. Early theories of 
cluster formation were the primordial- vorticity, pancake- 
sho ck, and tida l -torqu e theories (for a brief overview 
see iSueai fe Ivd Il995|) which had predicted different 
alignment configurations of cluster members for different 
birth models respectively. Unfortunately these models 
were at least incomplete due to the missing knowledge 
of the presence of dark matter in clusters and the 
assumption that clusters formed through fragmentation. 
Putting the focus on finding an orientation pattern for 
the whole cluster may be the main reason why many 
investigations in this field di d not yield sufficiently 
unambiguous resu lts l|Flinll2001l [Flin fc Godlowskil98fi 
ISueai fe Iv3ll995|) . 

More recent studies concentrated on the orientation 
of a number of brighter galaxies whose major axis 
shows a tendency to be aligned al ong the elongation 
of the parent cluster l | S trublel Il990) . For instance the 
Vir go l|W est fe Blakcslcc 20Q3) and also the Coma clus- 
ter l)Wesdl!997|l were found to exhibit this phenomenon 
which is furthe r discussed in sectio n [SJ 

Only recentlv lBrown et alJ l)2002|) published their work 
that tries to give an estimate on the reliability of grav- 
itational lensing studies which assume that anisotropies 
in the cllipticity distribution are due to a gravitation- 
ally induced shea r field distorting galaxy images (cf. 
iTvson et al J 1199 0). an important aspect briefly picked 
up in section [S] as well. 

In the present work, using extensive data on galaxy po- 
sition angles and inclinations of a large number of galax- 
ies in Coma (A1656), we are investigating the alignment 
of the majority of faint cluster members on small scales, 
even though only in a statistical way. Our method is to 
perform an individual statistical test for the surround- 
ings of each galaxy in the sample. We show that the 
found anisotropies are statistically significant to a high 
extent and superpose an image of average orientations of 



non-random regions on recent X-ray observations. Since 
an interesting coincidence emerges this approach seems 
to be profitable. We also try a tentative interpretation 
based on current models of galaxy cluster evolution. 

2. THE DATA 

A huge sample of 6724 galaxies within a regi on of 
ab out 2.6° x 2.6° centered on Co ma was compiled in l!983l 
bv iGodwin. Metcalfe, fc Peachl (GMP hereafter). This 
is still the most comprehensive catalogue of galaxies in 
the Coma cluster, comprising galaxies up to a bright- 
ness within the 26™5 contour limit of 626.5 = 21™0. The 
catalogue is considered complete up to 626.5 = 20™0. Po- 
sition angle and ellipticity for 4344 entries are included, 
however no morphological classification is given. 
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Fig. 1. — Distribution histogram of the 4344 position angles in 
the GMP sample. 

The sample as a whole shows a slight overabundance 
of galaxies with position angles around 90° (see Fig.^l 
with a roughly harmonic shape of period 180°. Also a 
periodicity of 45° may be inferred which was presum- 
ably caused by the plate scanning process but can be ne- 
glected for our purposes since it is smoothed out by the 
binning process. These glo bal properties hav e yet been 
known for quite some time l|Diorgovskilll988|) . However, 
a closer look at the respective orientations of the indi- 
vidual galaxies is worthwhile and provides new insights. 

3. THE METHOD OF ANALYSIS 

Instead of separating the cluster into certain sub- 
samples which requires us to make assumptions on where 
we expect to be alignment in the first place, we decided 
to use the environs of each galaxy as a sub-sample which 
eliminates such a selection bias. 

For each galaxy the angle of intersection between its 
major axis and the major axes of the 15 closest neigh- 
bor galaxies is calculated, whereby 15 was chosen as a 
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reasonable trade-off between small-number statistics and 
unwanted smoothing of the results. From these differen- 
tial angles a distribution histogram is created. If the 
distribution is not in agreement with a random distri- 
bution, i.e. it is not isotropic according to a x 2_ test, a 
best fit for the average orientation of this ensemble of 
adjacent galaxies is derived as described below. 

First, the 'projected spin vector' Pi of a galaxy with 
PA /3i and ellipticities q is defined as the vector in the 
observational plane that is orthogonal to the major axis 
of the galaxy image and has length Pi = sin (77). Here 
77 denotes inclination, thus cos (77) 
axial ratios -. 

a 

' COS (3i 

sin Pi 
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For the ensemble we define the 'average projected spin 
vector' S with PA a 



s = s( cosa ) 

V sin a I 

as the vector which satisfies the condition 
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The scalar product can be expanded to 

{S-P,) 2 = S 2 P 2 cos 2 {a ~p t ) (5) 

Thus P 2 can be viewed as the weighting factor of the sum 
which means that edge-on galaxies have a much stronger 
weight than (nearly) face-on galaxies for which the PA is 
more or less arbitrary. 
It follows that eq.0]has the solution: 

tan(2a) 



Si, 



£P?cos(2A) 



(6) 



If we identify x and y in eq.EJwe get the length of vector 
S by computing S 2 = S 2 + S 2 



= (^if cos(2ft)^ + (^if sin(2ft) 



(7) 



4. RESULTS 

The map seen in Fig.[5] shows the vector S plotted as 
black line at the position of each galaxy which is found 
to possess an anisotropic neighborhood with a level of 
confidence 7 X 2 = 0.99. It is interesting to note that no 
significant anisotropies are found around the central cD 
galaxies even though the galaxy density in this region 
is much higher than anywhere else on the map. More- 
over, clearly isolated groups of galaxies with anisotropic 
orientation can be found. 

Even though the algorithm favors such occurrence, the 
strength of these features is greater than could be ex- 
pected for a random sample. In order to give an esti- 
mate of the statistical significance of our results we per- 
formed a series of 1000 runs with artificially generated 
isotropic data. Apart from being impracticable due to 



Table 1. Definition of the samples 



sample 


brightness limit 


JV 


7 X 2 


1 


626.5 < 21 m 


4344 


0.99 


2 


&26. 5 < 20™ 


2354 


0.97 



the huge number of runs, a visual assessment of these 
artificial plots would also be subjective therefore we de- 
cided to apply a clustering algorithm to the distribution 
of S. We chose a modified version of K- means clustering 
which is simple to implement and fast enough for our 
purpose. The data were clustered in a three-dimensional 
space spanned by the x and y coordinates of S and as 
a third coordinate by its PA a. The original algorithm 
had to be modified to be cyclic in a. Also we used as a 
distance measure: 



Ax 2 + Ay 2 + Ctan 2 (Aa) 



(8) 



This definition makes sure that two S vectors which have 
a mutual PA difference of 90° are infinitely far apart as 
seen by the clustering algorithm. The scaling factor C 
can be viewed as weight of a. That is to say if C is large 
the clustering will be very sensitive to small differences 
in a and vice versa. 

The artificial galaxy samples were generated by assign- 
ing random PAs to each of the galaxies in the original 
sample. Moreover the inclinations were shuffled among 
the sample members which has the advantage of preserv- 
ing the inclination distribution. 

We performed two series of runs with 1000 artificial 
samples each. The first one included the whole set of 
4344 galaxies and a level of confidence for the x 2 -test 
of 7 X 2 = 0.99 was chosen while for the second run only 
those 2354 galaxies brighter than 6 2 6.5 < 20™0 were used 
which made it necessary to decrease 7 X 2 to a value of 
0.97 in order to get enough S vectors to work with (cf. 
Table P). 

The clusters produced by the K-means clustering had 
to be assessed in an objective way thus we introduced 
a quality value q for the features found in the artificial 
samples as well as in the original data: 

(a) Each feature that contains less than five members 
(i.e. S vectors) has q = 0. 

(b) For each remaining cluster q is defined as the norm 
of the sum over its constituent S vectors divided by 
its volume in the three-dimensional cluster space. 

Table shows the results obtained for samples 1 and 2. 
It allows to compare the number n and the quality q of 
the features in the data to the averaged values for these 
quantities calculated from the set of artificial samples. 
Here qA denotes the highest value for q found in the 
data whereas q~A is the average of the highest q values 
calculated from each individual artificial sample. 

U := 5>x)/1000 (9) 

The rightmost two columns in Table[21give the percent- 
age of artificial samples that showed smaller values for n 
and qA than the data. 

N( Xi < x) 



1000 



with x S {n, q A }, i e [1, 1000] (10) 
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Table 2. Properties of found anisotropies 



sample 




data 


artificial 


significance 


n 


QA 


n q A 


CV, CFq 


1 


51 


100000 


29 16900 


0.997 0.965 


2 


92 


42900 


71 10300 


0.957 0.951 



Table 3. Qualities of top three clusters of a sample 



sample 




data 




artificial 




er. 


QA 


QB 


QC 


QA QB 


qc 




1 


10 5 


95800 


4490 


16854 1054 


19 


1 


2 


42900 


11000 


10300 


10340 2391 


743 


0.999 



1000 - ^ . 



95800 - I 

'I} 

4490/ ' \ 

/ 




v 



ct„ := 



It becomes clear that only very few artificially 
generated samples showed more or better clustered 
anisotropies than the original data. 

In order to get a more sophisticated means of assess- 
ment we used the three best clusters in a sample de- 
noted A, B, and C in descending order of quality q. For 
the original data these would be the very strong cluster 
around NGC4911 (A) and the slightly less pronounced 
ones around NGC4839 (B,C). A detailed analysis re- 
veals that clusters A and B comprise 29 distinct galaxies 
each whereas 39 distinct galaxies cause the anisotropic 
region around cluster C. This can be seen in Fig.|3which 
also shows the spatial distribution of those galaxies. 

Table[3] gives the results from comparing the cluster 
qualities qA,qB,qc of the three top rated clusters in a 
sample. Here 

N ((& < q A ) A (q l B < q B ) A (gfc < q c )) 

1000 1 ' 

and i is the index of the artificial sample as above. None 
of the artificial samples could produce as many high qual- 
ity clusters of anisotropically distributed galaxies as the 
actual data. 

In order to understand their physical meaning it is 
desirable to verify the visually striking deviations from 
isotropy in Fig.|2] by comparing them with other mea- 
surements of anisotropy in the Coma cluster. An appro- 
priate means to this end can be the underlying image 
in Fig. |21 It shows the difference between the original 
X-ray intensity image of Coma and the expected inten- 
sity for a so called /3-model, which describes the distri- 
bution of hot intra-cluster gas in a completely relaxed 
cluster. This residual image was recently computed from 
a mosaic of XMM ob servations of Coma composed by 
INeumann et alJ l|2003|) who kindly put their data at our 
disposal prior to publication. According to them the pro- 
nounced feature in the SW around NGC4839 (white + 
symbol) is a subcluster currently falling on the main clus- 
ter. Due to ram pressure inflicted by the intra-cluster gas 
upon the subcluster, its gas is stripped off and leaves be- 
hind a track of hot gas. The same mechanism seems 
to have produced the slightly less strong residual around 
galaxies NGC4911 and NGC4921 located SE of the dom- 
inating cD galaxy pair in the center. 

In Fig. [3] indeed a striking coincidence is observed be- 
tween the regions of statistically significant galaxy align- 
ment and the features emerging in the residual im- 



x [arcsec] 

Fig. 2. — Map of vectors S for anisotropic galaxy subsamples 
(7 X 2=0.99). Dashed circles denote clusters identified by K-means 
algorithm, numbers mean quality q. The underlying image (ex- 
tension marked by light grey region) shows X-ray residuals gener- 
ated by subtracting a /3-modcl from the Coma cluster X-ray image. 
White + symbols denote massive galaxies. Origin at a=12 h 57™3; 
<5=+28°14.'4 (1950.0), north is up, east is left. 



age. Moreover it may also be worth noting that the 
S vectors around NGC4839 seem to point in a direc- 
tion tangential to its perimeter while in the case of the 
NGC4911/NGC4921 pair the anisotropies seem to coin- 
cide with the connecting line between the two galaxies. 

5. INTERPRETATION 

A possible explanation for the intriguing match found 
in Fig.[5]and Fig-El might be given by a simple model of 
the dynamical evolution of interacting galaxies falling on 
Coma along a filam ent, a process seemingly very common 
for galaxy clusters llKim et alJl2002T: iNovikov et al.|[T999t 
Plionis & Basilakos 2002). What we called the Laminar 
Flow Model (LFM) is outlined in Fig.|J 

(a) Galaxies are streaming into the main cluster from 
outside. The velocity vectors of two neighboring 
galaxies can be decomposed into components par- 
allel («||) and perpendicular (v±) to the velocity 
vector v of the center of mass (CoM). 

(b) If viewed in a co-moving reference frame with ori- 
gin in the CoM, the two galaxies orbit each other 
prior to merging. Their angular momentum vector 
L ~ r X Vj_ lies in a plane orthogonal to v± by 
definition. Thus the projection of L upon the ob- 
servational plane is statistically more likely to be 
aligned parallel to the infall direction v. 

(c) Finally, after the merging has taken place, the re- 
sulting galaxy tends to possess a spin vector paral- 
lel to the previous angular momentum L. 

Basically the process can be viewed as funneling of galax- 
ies along the filaments. Perpendicular to the infall di- 
rection matter is coalescing and getting denser whereas 
parallel to its motion vector v it is stretched due to tidal 
forces. Thus the probability for a merger is higher if the 
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Fig. 3. — Enlargements from Fig.l2l showing the best three clusters A,B,C and their surroundings. Black dots denote those galaxies 
which comprise the neighborhood of the S vectors in the cluster. 




Fig. 4. — The Laminar Flow Model: (a) Stream of infalling 
galaxies - (b) Orbiting galaxies - (c) Merger product 



component of the radius vector r parallel to v is small. 
From the definition of L and Fig.0{b) it follows that the 
observed spin alignment can be explained by the LFM. 
Moreover it is i n acco rdance with the view favored by 
IMaller fe Dekell <|2002f) which holds that the spin vector 
of galaxies is chiefly determined by the last few merging 
events. 

The os tensible contradict i on wi th the findings of, for 
instance, iWest fc Blakesleei (|200f)D that the major axes 
of the brightest cluster galaxies are aligned with the fila- 
ments, can be resolved by considering the different mech- 
anism that produces the bright central cD galaxies. They 
are the final destination of the infalling matter whose lin- 
ear momentum is thus converted into angular momentum 



by the ultimate merger with one of those cosmic canni- 
bals. 

Finally the issue of weak lensing should not remain 
unmcntioncd. On one hand a detected orientation 
anisotropy in a sample of galaxies may be caused by grav- 
itational lensing, on the other hand if the anisotropics are 
real they may give wrong lensing indications if not al- 
lowed for. These considerations are beyond the scope of 
this paper but further investigation may be worthwhile. 



6. CONCLUSION 

In this work we find that the overall isotropic appear- 
ance of galaxy orientations in Coma can not be main- 
tained when looking at galaxy ensembles on the smaller 
scales of subclusters. These fluctuations are not caused 
by statistical variations but there exists strong evidence 
that they are the result of anisotropic merging of subclus- 
ter members which fall on the main cluster presumably 
along filaments extending between large clusters. There- 
fore the identification of regions in which galaxies show 
statistically significant deviations of their spin vectors 
from isotropy may help to trace such filaments as well as 
provide a new and useful tool to investigate the evolution 
of galaxy clusters. 

Acknowledgments: The mosaic of X-ray observations of 
the Coma cluster was provided by courtesy of D. Neu- 
mann fcf. lNeumann et al.ll200^ . 
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